ON THE CLARK-a MODEL OF TURBULENCE: 
GLOBAL REGULARITY AND LONG TIME DYNAMICS 



CHONGSHENG CAO, DARRYL D. HOLM, AND EDRISS S. TITI 

Abstract. In this paper wc study a well-known three— dimensional turbulence model, the filtered Clark 
model, or Clark— a model [7] . This is Large Eddy Simulation (LES) tensor-diffusivity model of turbulent 
flows with an additional spatial filter of width a. We show the global well-posedness of this model with 
constant Navier-Stokes (eddy) viscosity. Moreover, we establish the existence of a finite dimensional 
global attractor for this dissipative evolution system, and we provide an anaytical estimate for its 
fractal and Hausdorff dimensions. Our estimate is proportional to (L/l^) 3 , where L is the integral 
spatial scale and l d is the viscous dissipation length scale. This explicit bound is consistent with the 
physical estimate for the number of degrees of freedom based on heuristic arguments. Using semi- 
rigorous physical arguments we show that the inertial range of the energy spectrum for the Clark-o 
model has the usual A; -5 / 3 Kolmogorov power law for wave numbers ka -C 1 and fc -3 decay power law 
for ka 2> 1. This is evidence that the Clark— a model parameterizes efficiently the large wave numbers 
within the inertial range, ka 2> 1, so that they contain much less translational kinetic energy than their 
counterparts in the Navier-Stokes equations. 



1. Introduction 

This paper is devoted to the mathematical analysis of a Large Eddy Simulation (LES) model of 
turbulent flows with explicit filtering. We consider the "filtered Clark model," aka "Clark— a model" [7]. 
This is a nonlinear LES model of turbulence consisting of the "tensor-diffusivity model" of Leonard [19], 
filtered by inversion of the Helmholtz operator with width a. As pointed out in Winckelmans ct al. [24], 
this model is generic: For all regular symmetric filters that have a nonzero second moment, this form 
is found as the first term of the reconstruction series for the filtered-scale stress. Thus, the Clark— a 
model we consider is generic and corresponds to significant, yet perhaps incomplete, reconstruction of the 
filtered-scale stress. Such an explicitly filtered tensor-diffusivity model was used by Vreman et al. [22, 23] 
with good success, in combination with a dynamic Smagorinsky term as a "nonlinear mixed model." In 
their work, the Smagorinsky term was posed to model the truncation effects of subgrid-scalc stress and 
incomplete reconstruction of the filtered-scale stress. Vreman et al. [22, 23] referred to their model as 
"a dynamic version of the mixed Clark model," while Winckelmans ct al. [24] wrote that they preferred 
the term "tensor-diffusivity model" for such models. The results of [24] supported the view that the 
explicitly filtered tensor-diffusivity (or, Clark) model may well suffice for practical reconstruction of the 
filtered-scale stress, without introducing a Smagorinsky term. 

Section 2 is devoted to the precise definition of the Clark— a model and to discussions of the properties 
of its nonlinear terms. Section 3 establishes the main results of global (in time) regularity for the 
Clark— a model. That is, global existence, uniqueness and continuous dependence of the initial data of 
the solutions to this 3 — D model of turbulence. Section 4 shows that the "Clark— a turbulence" model 
possesses a global attractor and gives an upper bound on its Hausdorff and fractal dimensions. This 
upper bound is discovered to scale as (L/ld) 3 - the cube of the integral scale L divided by the dissipation 
length Id = (i^/e) 1 / 4 . Finally, the translational kinetic energy spectrum E{k) for the Helmholtz-filtered 
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Clark— a model is also shown in section 4 to pass from fc~ 5 / 3 for smaller wave numbers (ka <C 1), to fc~ 3 
for larger wave numbers (ka 3> 1). It is worth mentioning that the results presented here are similar 
to those already established for the Navier-Stokes-a (NS-a) model [11], [12] (also known as the three 
dimensional viscous Camassa-Holm equations, or the Lagrangian averaged NS-a model). Furthermore, 
when viewed as a closure model of turbulence for the Reynolds averaged equations in turbulent channels 
and pipes the Clark-a model gives exactly the same reduced equations in channels and pipes as the 
NS-a model. Comparing the solutions of those reduced equations in pipes and channels to the empirical 
data has already given excellent agreement (sec, [2]— [5].) Therefore, in this regard, one can assert that 
the Clark-a model is as successful as the NS-a sub-grid model of turbulence. 

2. Preliminaries and Notations 
Let fl = [0, 2ttL} 3 . The so-called Clark— a turbulence model (cf. [7]) of viscous incompressible flows 



in a domain Q subject to the periodic boundary condition reads: 

d t u - vAu + (u • V)u + Vp - V • (H _1 (a 2 Vw • Vw T )) = /, (1) 

V-u = 0, (2) 

u(0) = u , (3) 



where, u represents the unknown "averaged/filtered" fluid velocity vector, and p is the unknown "av- 
eraged/filtered" pressure scalar; v > is the constant kinematic (eddy) viscosity, a is a length scale 
parameter which represents the width of the filter. The body forcing term, /, which is time independent, 
and the initial velocity, uo, are given. The operator Ti. is the Helmholtz operator, defined as 

TLu = u — a 2 Ait, subject to periodic boundary condition, 

and the tensor (Vu • Vu T ) is given by 

3 

(Vu • Vji T ) tJ = Vu t ■ Vu 3 = ^2(d Xk Ui) (d Xk Uj). 

k=l 

From now on, the Einstein's summation convention will be used, e.g., 

3 

{d Xk Ui)(d Xk Uj) = ^(d^Ui^d^uj). 

k=l 

Let v = H.u = u — a 2 Au. Then, the above system is equivalent to 



d t v - uAv + (u ■ V)v + (v ■ V)u - (u ■ V)u - a 2 V • (Vit • Vit T ) + X7q = g, (4) 

V-u = V-w = 0, (5) 

v(0) = v = Hu 0l (6) 
where, q = Tip and g = Ttf. 



2.1. Notations. Let £ p (f2) and H m (n) denote the usual L p Lebesgue spaces and Sobolev spaces, re- 
spectively (cf. [1]). We will denote by | • | the L 2 — norm, and (•, •) the L 2 — inner product. Let T be the 
function space which consists of all the vector trigonometric polynomials. We set 

V = {0eJ ?r :V-(/) = O & J <t> = Q}- 

The spaces H, V\ and Vi will denote the closures of V in L 2 .H 1 and H 2 respectively. Let P a : L 2 — > H, 
be the orthogonal projection, and let A = —Pa- A be the Stokes operator subject to the periodic boundary 
conditions. It is well known that A^ 1 is a self-adjoint positive compact operator in H and V(A) = V2 
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(cf. [9], [15], [17], [20]). Let < Ai = L 2 < \ 2 < ■ ■ ■ be the eigenvalues of A, repeated according to 
their multiplicities. It is well-known that a constant C\ > exists such that 

^ < ^ < C l3 2 ' 3 , (7) 

for j = 1,2, • • • , (see, e.g., [9]). 
For w\,w 2 G V, we define 

B(w u w 2 ) = P a ({ Wl -V)w 2 ). 
The bilinear form B has the following properties (cf . [9] , [20] ) . 

Proposition 1. The bilinear form B can be extended as a continuous map B : Vi X Vi — > V{, where V{ 
is the dual space of V\ . Furthermore, 

(B(w 1 ,w 2 ),w 3 ) v , = - (B(w 1 ,w 3 ),w 2 ) v , 

for every w\,w 2 ,w 3 G V\. 

Proposition 2. Let u £ H 3 R V\, v = Hu = u — a 2 Ait. Then, the bilinear operators B(u, v), B(v, u) and 
3 

^^B(d Xj u,d Xj u) are well defined with values in V{. 
i=i 

The following properties can be established by integration by parts. 

Proposition 3. Let u G H 3 PI Vi, v = Hu = u — a 2 Aw. Then, 

(i) (B(u,u),u) v , =0, 

(ii) (B(v,u),u) v , = 0, 



(iii) (B(u,v),u) v , - a 2 ^2 (B{d X] u,d X] u),u) v , = 0. 



3=1 

As for the case of the Navier-Stokes equations (cf. [9], [20]) by applying the operator P a , we get the 
following system which is equivalent to the system (l)-(3) 

du 

— + vAu + B(u, u) + a 2 7i- 1 B{d X] u, 8 Xj u) = f, (8) 
u(0) = u , (9) 
and the following system which is equivalent to the system (4)-(6) 

^ + vAv + B(u, v) + B(v, u) - B{u, u) - a 2 B(d Xj u, d Xj u) = g, (10) 
v(0) = v . (11) 

Definition 4. Let T be any fixed positive number and f G Vi . We call u a regular solution of (8)-(9) 
as an equation in V± on [0, T] if u satisfies the system (8)-(9) and 



ue C w {[0,T},V 2 ) n L 2 ([Q,T], H 3 (n)), and 

^eL 2 ([o,n Vl ). 

Here C w ([0, T], V 2 ) is the functional space of all weakly continuous functions from [0, T] to V 2 . 
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For convenience, we recall the following Sobolev and Ladyzhenskaya's inequalities in M 3 

M^KCoMli'M 1 ^, (12) 

\\u\W<C Q \\ut^\\u\\%t, (13) 

HU»<c |H|fli, (u) 

for every u € V\ , and the Agmon's incqulity in R 3 

NU~<Coh||^ 2 |H|^ 2 , (is) 

for every u S Vi- Here Cq is a universal constant. Also, notice that 

\Hu\ 2 < (Ar 1 + a 2 ) (|Vu| 2 + a 2 \Au\ 2 ). (16) 
For this model the quantity 

(«,«) = |u| 2 + a 2 |Vu| 2 

will be called the energy of the system, while the quantity [u, u) — \u\ 2 will be called the translational 
energy of the system. 

3. Existence and Uniqueness 

3.1. H 1 estimate. The estimates presented here are formal. They can be made rigorous by establishing 
them first for a Galcrkin approximating scheme and then passing to the limit by using Aubin compactness 
theorem (see, e.g., [9], [20] and [21]). By taking the V{ action of equation (8) with v = Hu and using 
Proposition 3 and Lemma 1.2 in Chapter III of [20], one finds the energy balance satisfies: 

1 d(\n\ 2 4- rr 2 Wii\ 2 ) 

3 U ' ° ' ' ) + "(\Vu\ 2 + a 2 \Au\ 2 ) + (B(u,u),Hu) v , - a 2 (B{d Xj u,d Xj u),u) v , = (f,Hu) v , . (17) 
By using part (iii) of Proposition 3 we obtain 

1 d(\u\ 2 + a 2 |Vu| 2 ) 9| . l9 . ,,„, \ 

2 dt + K|Vu| 2 + a 2 \Au\ 2 ) = (f,Hu) (18) 

< |/| |u| + a 2 |V/||Vu| (19) 

< ^(l/l 2 + a 2 |V/| 2 ) + ^(\u\ 2 + a 2 \Vu\ 2 ) (20) 

< ^d/l 2 + « 2 |V/| 2 ) + ^(|Vu| 2 + a 2 \Au\ 2 ). (21) 



As a result, we reach 



d(|u| 2 +a 2 |Vu| 2 ) „ 2|A |2N |/| 2 + a 2 |V/| 2 

-j—! — + ^(|Vu| 2 + a 2 |Au| 2 ) < m 1 - 1 . (22) 

Thanks to the Gronwall inequality, we get 



u 



2 1 a 2 \Vu\ 2 < ATi(Q!,i/,t), (23) 



where 

#i(a, i/, t) = e-" Alt (| U o| 2 + a 2 |Vw | 2 ) + - ' ' J ' ■ (24) 



Jl 2 + « 2 |v/l 2 

j/ 2 A? 



Moreover, by (22), for each r > 0, we have 

H f |2 i „,2|V7 f|2 



-/ r (|V M | 2 + a 2 |A^| 2 ) < +Q ' |V/I )r +g 1 (a, I /,t). (25) 
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Consequently, 



2 IA ,2^ (|/| 2 +« 2 |V/ | 2 )^ , 1 



(|Vu| 2 +a 2 |Au| 2 ) < 



+ -(K| 2 + a 2 |Vw P 



(26) 



3.2. H 2 estimate. Again, the estimates presented here are formal. As above, they can be made rigorous 
by establishing them first for a Galerkin approximating scheme and then passing to the limit. By taking 
the V{ action of equation (10) with v = Hu, one finds 

\^dt + ^ VV ' 2 + ( B ( u > v )+ B ( v > u ) ~ B ( u ' u ) ~ a 2 B(d Xj u,d X] u),v) v , = {g,v) v , . (27) 



Notice that 



(B(u,v) + B(v,u) - B(u,u) - a 2 B(d X] u, d Xj u),v) v , 

= (B(v,u) - B(u,u) - a 2 B(d Xj u, d Xj u) ,v) v , 

< |V«| ||«||i4 + \v\ \\u\\ L s\\u\\ Le + a 2 \Vv\ \\Vu\\ 2 Li 

< C 2 \Wu\ {v^lW^ + C^vl IuI^IVuI^ + CqVIVuI \Vu\^ 2 \Au\ 3 / 2 . 



Also, 



(9,v)vi = W,v) v , = (f,v) v , - a 2 (Af,v) v , 
< l/l K'I+« 2 |V/| |V«|. 
By Cauchy-Schwarz and Young's inequalities we reach 
d\v\ 2 v 



dt 4 



|V«| 



< l/l 2 + M 2 + « 4 |v/| 2 + 2C ^f^ + sc 2 \ v \ WW + 8C ^ ' Am ' : 



< |/| 2 + a 4 |V/' 2 



2Cg|Vu| 4 SCgAyVl^lVul 1 / 2 8Cff|VM| |Au| 

1 "T a + o 



By (23), (26), and Gronwall inequality we get 

KOI 2 < K 2 (a, u, t) (\v\ 2 (s) + t(\f\ 2 + a 4 |V/| 2 )) , 

where 

72 \ 3/4/ 



K 2 (a, v, t) = exp 



i | 2C$(K 1 (a,v,0)) 2 t SOn^^ | 8C 4 |/| 



u 3 a 2 X 1 



1/2 



+ 



8QK| 2 +a 2 |Vu | 



In particular, 
Moreover, 



K*)| 2 < K 2 {a,v,t) {\v \ 2 + t(\f\ 2 + a 4 |V/| 2 )) 



i/ 

4 7o 



|V V | 2 < (|/| 2 + a 4 \ V/| 2 ) i + A 2 (a, ^ t) (M 2 + i(|/| 2 + « 4 |V/| 2 )) 
These analytical estimates lead to the following. 



(28) 



(29) 
(30) 
(31) 
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Theorem 5. Let uq G V\ and f £ Vi. Then there is a unique regular solution to system (8)-(9) for 
all t > 0. Furthermore, this solution depends continuously on the initial data in the sense that it will 
specified in the proof below. 

Proof. One can establish the existence of a regular solution by applying the standard Galerkin approx- 
imation procedure together with the a priori estimates (23), (26), (30) and (31) (see, for example, [9], 
[15], [17], [20]). Here, we shall only show the uniqueness of regular solution to the system (8)-(9) and the 
continuous dependence of solutions on the initial data. Suppose that u\ and u 2 are two regular solutions 
to the system (8)-(9), and let w — u\ — u 2 . Then the difference w satisfies 

^ + vAw + +B(w,u 1 ) + B(u 2 ,w) + a 2 n- 1 B{d X] w,d X] u l ) + a 2 n- l B{d X3 u 2l d X] w) = 0, (32) 
w(0) = 0. (33) 

By taking the V{ action of equation (32) with Tiw and applying Lemma 1.2 in Chapter III of [20], one 
finds 

1 d(UnV -I- n, 2 \XIii,\ 2 \ 

v{\Vw\ 2 + a 2 \/\w\ 2 ) 



2 dt 

= -(B(w,ux) +B(u 2 ,w) + a 2 U~ 1 B(d Xj w,d Xj u 1 ) + a 2 H~ 1 B(d Xi u 27 d X] w),Hw) v , 

< |Wifl|(|M|L-|Vui| + IMl~|VH) +a 2 (|Viii + |Vu 2 |)||Vu;||i< 

< C Q \Hw\{\Vw\ 1 ' 2 \^w\ 1 l 2 \Vu l \ + |Au 2 ||VH) +a 2 C {\V Ul \ + |V?i 2 |)|V«;| 1/2 |AH 3/2 - 
By Cauchy-Schwarz inequality, we obtain 

1 d(\w\ 2 + a 2 \Vw\ 2 ) ^ ,. . ,„ ,„ 

2^ — dt — - 2 ' + 1 11 + 1 1 1 ' 

where C is a constant. Thanks to the Gronwall inequality, we get 

\w{t)\ 2 + a 2 \Vw(t)\ 2 < exp jf (|A M2 | 2 + |V Wl | 4 + |Vu 2 | 4 )^) (\w(0)\ 2 + a 2 |Vu;(0)| 2 ) . 

Thanks to (23), (26), (30) and (31), the above implies continuous dependece on initial data. In particular, 

«;(*) = when to(0) = 0. 
Therefore, the regular solution is unique. □ 

Remark 1. In order to be able to estimate the dimension of the global attractor it is required that the 
solution is differentiable, in the appropriate norms, with respect to the initial data. Following similar 
energy methods to the ones introduced above one can establish the required differentiability (see, e.g., 
[8], [9] and [21]). 

4. Global Attractors and Energy Spectra 

In this section we show the existence of the global attractor. Moreover, we provide an upper bound 
to its fractal and Hausdorff dimension. Also, we consider the energy spectrum for the Clark— a model. 

4.1. Global Attractor. Denote by u(t) — S(t)uo the solution of the system (8)-(9) with initial data 
uq. As a result of Theorem 5, one can show that 

u(t) = S(t)u a e V x for all u G Vi,t > 0, 

and 

u(t) = S(t)u G V 2 for all u G V 2 ,t > 0. 
Since, in this section, we only consider the long time behavior of solutions of the system (8)-(9), by (26), 
(28) and Theorem 5, we conclude that u(t) G Lg c ((0, S},V 2 ) for every u(0) G V\ and any S > 0. 
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Theorem 6. Suppose that f G V\. Then, there is a compact global attractor A C V\ for the system 
(8)-(9). Moreover, A has finite Hausdorff and fractal dimensions. 

Proof. First, let us show that there is an absorbing ball in V\ and V2. Let u be the solution of the system 
(8)-(9) with initial data it 6 V\ and u | 2 + a 2 |Vu | 2 < p. By (22), we have 

Kt)| 2 +a 2 |V^)| 2 < e- v ^\\utf + a 2 |Vu | 2 ) + lfl " + f^ f ^ 

" 9 + ^ 2 A 2 

As a result of the above, when t is large enough such that 

P " ^ 2 A 2 

we have 

| U (t)| 2 + a 2 |V U (t)| 2 <R a {a,v, f), (34) 



where 



In particular, 



,I/I 2 + « 2 |V/I 2 
v 2 \\ 



R a (a,v,f) = 2 rT "2 . (35) 



limsup(Hi)| 2 +a 2 |V U (£)| 2 )<2l /|2 



t-00 ^ 2 Af 
Therefore, system (8)-(9) has an absorbing ball B in V\ with radius R a (a, v, /). 

Next, we show that there is an absorbing ball in V2. First, by (34) and (25) we have 

fW + «»|A„|») < „ fW + Q =|A„|») < l/l' + + ">*■>, (30) 

Jt Jt ^ A x 

By applying the uniform Gronwall inequality (cf. for example, [21], p. 89) and (36) and (26), we obtain, 
when t is large enough, 

\Hu(t)\ 2 = \v(t)\ 2 <R v (r, a, u,f), (37) 

where 



R v (r,a,is,f) = exp 



'r 1 / 2 



««,">/)) + ^7*^ ^ ) J 



'(|/| a + a 2 |V/| 2 )(2 + ri/Ai) , r lr 



i^ 3 A 2 r 



-IV/N, (38) 



and r > is fixed. Therefore, we conclude there is an absorbing ball B in V2 with radius R v (r, a, v, /). 
Thanks to Rcllich Lemma (sec, e.g., [1]), the operator Sit) is a compact operator from V\ to itself. 
Following the standard procedure (cf., for example, [8], [9], [10], [18], [21] for details), one can prove that 
there is a global attractor 

A = n s>0 (U t>s S(t)B)cV 2 . 

Moreover, A is compact in V\. 
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By using the estimates (34) and (36), applying Lemma 4 and Lemma 5 in [11], and following the proof 
of Theorem 6 in [11] line by line, we obtain the Hausdorff and fractal dimensions of the attractor A 

, (M< , (M<r [ Q/I 2 + " 2 |V/I 2 ) 2/3 (l/l 2 + " 2 |V/l 2 ) 3/4 ) 
d H (A) < dpi A) < C max <! ^ri > tft } ■ 

Moreover, following [11], we define the mean rate of dissipation of energy, 

e = suplimsup- / (|| V M (t)||| 2 + a 2 \\Au{t)\\ 2 L2 ), (39) 

u £A T^oo J Jo 

where u(t) is the solution corresponding to the initial datum uq. By analogy, we define the Kolmogorov 
dissipation length scale corresponding to the system (l)-(3), 

u - 

Hence, again, following the proof of Theorem 7 in [11], we have 

L\ 3/4 /L 



d H (A) < d F (A) < C - , 

We refer the reader to [11] for more details. □ 

5. Spectral Scaling 
Following [12] (see also [13] and [14]), we set 

1 





(2ttL) 3 


/ u[x)e 
Jq 


~ tK - x dx, 


Vk 




/ v(x)e~ 

JQ 


- lk ' x dx, 


Uk 


= E 

fe<|j|<2fc 


Uje ij ' x , 




Vk 


= E 

fc<|j|<2fc 






< 


= E^> 

j<k 


= 


E v j 




= E u ^ 

2k<j 




= E V 3 

2k<j 



The energy balance equation for the Clark— a model for an eddy of the size k 1 is 

\ li^ Uk,Vk ^ + ^(~ A " fc ' Wfc ) = Tk _T 2fe, (40) 
where the right hand side is the energy flux and 

T k = - {B(u<,u<),v k )~a 2 (B(Vu<,Vu<),u k ) + 

+ (B(u k + u>,u k + u£), v<) + a 2 (B(Vu k + Vu>,Vu k + Vu>), u<) . 
Taking an ensemble average (long time average) of (40) we get 

(v(-Au k ,v k )) = (T k ) - (T 2k ). (41) 
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Let 

£ Q (fc) = (i + c* 2 |fc| 2 ) J2 M 2 - 

\j\=k 

Then (41) can be written as 

,-2k 

vk 3 E a (k) ~ / k 2 E a (k) dk ~ (T fe ) - (T 2k ) . 

Jk 

For wave numbers k within the inertial range it is assumed that there is no energy disspiation, hence we 
have vk 3 E a {k) <C (Tfc), and we also have (Tfc) w (T2k), because there is no leakage of energy. One can 
consider three possible scales for the average velocity on an eddy of length size fc _1 . Namely, 

^ 0) = (^3 Jjv k \ 2 d x y 2 ~ (k(l + a 2 k*)E a (k)) 1/2 , 
Ujp = ^ J^u k -v k dx\ ~ (kE a (k)) 1/2 , 



1 f , l2j \ 1/2 / fcS Q (fc) \ 1/2 



ul 2) = (4-J W?d X 



L 3 J a [ I \l + a 2 k\ 

The corresponding turnover time T k for an eddy of the size k~ l will be (cf. [16]) 

1 (lWfc 2 )^- 1 '/ 2 



k kUl fc 3 / 2 (£«(*;)) 1/2 ' 
The energy dissipation rate e (39) is 



0,1,2. 



1 ? k Ej m *•"«.(»»"■ 



r. 



(l + a 2 ^ 2 )^- 1 )/ 2 ' 



e^^l + a 2 ^ 2 )^- 1 )/ 3 



As a result, we have 

Therefore, the translational kinetic energy spectrum E(k) = X)|j|=fc "j'| 2 ^ s given by 

f e 2/3 

Fm= ^q(fc) F73' when to «1, 

— ^ — when fca » 1 . 

I a 2(4-n)/3^.(13-2n)/3 

We refer the reader to [6] and [12] for more details and discussions of the implication of the energy 
spectrum fc-( 13 ^ 2 ™)/ 3 for the larger wavenumbers (1 <C ka). 

(n) 

Remark 2. It is worth stressing that the filtered Clark model does not specify which one of the Up, 
for n = 0, 1, 2, is the right average velocity for an eddy of the size fc _1 . Consequently, it is not evident 
what would be the correct energy spectra in the sub-range, ka ;§> 1, of the inertial range. This will 

(2) 

be a subject of future research. However, our earlier research suggests the choice Ul , for which the 
translational energy spectrum for the filtered Clark model has the usual k~ 5 / 3 Kolmogorov power law 
for wave numbers ka -C 1 and shows fc -3 decay power law for /;« » 1. 
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